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Abstract: 

A physical nonlinear dynamical model of a laser diode is considered. We propose a feed-forward 
control scheme based on differential flatness for the design of input-current modulations to 
compensate diode distortions. The goal is to transform without distortion a radio-frequency 
current modulation into a light modulation leaving the laser-diode and entering an optic fiber. 
We prove that standard physical dynamical models based on dynamical electron and photons 
balance are flat systems when the current is considered as control input, the flat output being the 
photon number (proportional to the light power) . We prove that input-current is an affine map 
of the flat output, its logarithm and their time-derivatives up to order two. When the flat output 
is an almost harmonic signal with slowly varying amplitude and phase, these derivatives admit 
precise analytic approximations. It is then possible to design simple analogue electronic circuits 
to code approximations of the nonlinear computations required by our flatness-based approach. 
Simulations with the parameters of a commercial diode illustrate the practical interest of this 
pre-compensation scheme and its robustness versus modelling and analogue implementation 
errors. 
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1. INTRODUCTION 

Radio-over-Fiber is a very attractive technique for wireless 
access network infrastructure, because it can transmit mi- 
crowaves and millimeter-waves through optical fibers for a 
long distance. Therefore, the distribution of radio signals 
over optical fiber, which is of great interest for many 
applications such as broad-band wireless access networks, 
sensor networks, radar and satellite communications, has 
been intensively studied. And there have been rapid ad- 
vances in the techniques to generate and transport radio 
signals over optical fiber in recent years (see, e.g., Ai- 
Raweshidy and Komaki [2002]). When the wireless channel 
is in series with the optical link, nonlinear distortion due 
to the electrical/optical conversion is the biggest concern 
(see, e.g., Maury et al. [1997]). This makes nonlinearity 
compensation an attractive solution to improve link per- 
formance. 

The focus of this work is to investigate the standard phys- 
ical model used for describing commercial laser diodes. 
This model is based on rate equation for the electron 
and photon populations. This allows us to propose a pre- 
compensation based on differential flatness for the design 
of the input current in order to cancel diode distortions. 



Simulations corresponding to a commercial laser diode 
illustrate the pre-compensation scheme. When the distor- 
tion due to the optic fiber is negligible, our flatness-based 
pre-compensation scheme can be interesting to increase 
the bandwidth and transmission rate. 

This paper is organized as follows. In section 2 we recall 
the standard diode model and show that it is structurally 
flat, with the photon number being the flat output y. We 
show via elementary computations that the input current 
is a linear combination of y, 4r<y, 4pV, Jj(logy) and 

^2 (logy) with coefficients depending via explicit formula 
on the physical parameters (equation (3)). This relation is 
the starting point of the flatness based pre-compensation 
scheme addressed in section 3. In Section 4, we propose 
natural approximations when y is quasi-harmonic with 
slowly varying amplitude and phase. These approxima- 
tions simplify considerably the computation burden in 
such a way that a nonlinear analogue circuit can be eas- 
ily design to realize our flatness-based pre-compensation 
scheme. In conclusion, we propose some hint to deal with 
phase influence that becomes important when fiber distor- 
tions cannot be negligible. 
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2. THE FLAT PHYSICAL NONLINEAR MODEL 

A standard nonlinear model relating the input-current / 
to the power of the emitted light P is recalled in Lee et al. 
[2003] where this model is used to analyze the nonlinear 
"transfer" between input-current modulation / and the 
resulting output P) (see also Agrawal [1997] for a simple 
exposure on such basic models): 



P = R(P, X) 



d_ 

dt 
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where 



• t p and r n are photon and carrier lifetimes; 

• X is the normalized carrier density; 

• I t h is the threshold current and a is the constant 
related to other physical parameters, a = FI t hT n /T p . 

• R(P, X) is the net rate of stimulated emission. Fol- 
lowing Lee et al. [2003] we use 



R(P,X) 



' BT n I th {X - 1) + ^ 

1 + FBt p t c P 



where B, t„,t c are physical parameters . 

Since the input-current appears only in the second differ- 
ential equation of (1) (the electron balance equation), the 
first state, P, is a flat output. More precisely, this means 
that, if we consider (1) as a control system with input / and 
output y = P, we have a nonlinear input/output system 
such that its inverse admits no dynamics. This system is 
differentially flat (see Fliess et al. [1992, 1995], Martin 
et al. [2003], Sira- Ramirez and Agarwal [2004]) with P 
as flat output: the input / is a nonlinear function of P, 

^f, jpP- Notice that this property is independent of the 
precise form of R(P, X) when it depends effectively on X. 
This is always the case for physical reasons. 

Let us detail now the computations with the precise form 
of R given here above. The first equation of (1) yields 

v i P d f P \ d , ( P 

where Pi — f ^ 7 T n and t\ = B ^ lth are constant parameters. 
Thus we can derive this expression versus t to get |Xas 
a function of P, -j^P and J^-P: 



lx = -(- 

dt dt \P[ 



d 2 [P\ d\ ( P 



A linear combination of the two equations of (1) gives / 
explicitly: 

7— = X + r n X + -^-(P + P/t p ). 

J-th T n I t h 

Thus we have the explicit formula relating / to the 
derivatives of P: 

P\ ( i , T„T„\ d (P 
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With the normalized input/ouput variables 



dt \Pi 
P 



u = 
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Pi 



(2) 



we have the following simple but nonlinear input/output 
relationships that is equivalent to the physical model (1): 



u= 1+ 1+ 



' n ' p 



dV 



d 2 ( d n d 2 n \ , oX 

+ r n r p —y + T l ^-(\o g y)+r n — (lo g y)j (3) 

where the four parameters r„, t p , t c and r; = Br i th 
are positive time-scales. Remember that the normalized 
electron density A is a combination of y and ^y: 

d d_ 



X = 1 + y + t p — y + n — (logy). 



3. FLATNESS-BASED PRE-COMPENSATION 

We can now use (3) for deriving feed-forward strategy in 
order to compensate the non-linear distortion due to the 
diode. Assume that our goal is to transfer the amplitude 
and phase modulated signal e(t) cos(ivt + <fi(t)) from the 
current to the light leaving the diode. Here u> is the carriage 
pulsation, typically around a few Ghz for radio frequencies, 
and e(t) and <j>(t) are slowly varying amplitude and phase 
where the bits are encoded. This means that |^e| <C we 
and *C ui. Without lost a generality, we can assume 

that e < 1. Assume finally that this modulated signal has 
to be converted in photons y > according to y = y(l + 
e(t) cos(uit + (j>(t))) where y is some positive constant. This 
means that our goal is to find the current u leading to 
such y. If one knows the diode parameters we can use (3) 
to compute u. 

We have tested in simulation this simple pre-compcnsation 
scheme for the commercial diode considered in Lee et al. 
[2003]. For this diode, we have 

T n = 179 ps, Tp — 4.33ps, t c = 3.18 ps, ti = 1.81 ps. 

We choose a modulation frequency ui of 10 Ghz (w = 
2tiT0 10 s- 1 ), <j) = and 



e(t) = 1(1- (■( is 



(4) 



The average normalized light power is y = 0.0175 corre- 
sponds physically to P = 1.4 mW since Pi — 78.5 mW. 
The simulations of figure 1 illustrate the interest of taking 
into account the derivative terms in (3) to modulate the 
input current u. 

4. ANALOGUE IMPLEMENTATION ISSUES 

For such high carrier frequency u <~ 10 GHz, it is 
difficult and almost impossible to use real-time numerical 
computations. Let us now propose two approximations 
such that the above pre-compensation scheme can be 
applied in real-time via a specific electronic circuit. 

The first approximation is relative to the computation 
of derivative of y. Since we are looking for u such that 
y = y(l + e{t) cos(u;t + <fi(t)) with e and <fi slowly varying 
quantities, we have 
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Fig. 1. Input current modulations with e(t) given by (4): 
the solid curve corresponds to the response of the 
diode when u is related to y = y(l + e(i) cos(wi)) 
via (3); dashed line corresponds to the response of 



dt 



the diode when u = 1 
(static model). 

y ~ — yewsin(W + cf>), 



1 



dt 



y(l + e(t) cos(wt)) 



: y « — yeus 2 cos(cjt + < 



(5) 

The second approximation is relative to the nonlinearity 
attached to the log-terms in (3). Now we will use £< 1 
and propose an approximation up to second order terms: 

-^-(logy) ~ —eLusmtiut + 0)(1 — ecositut + 4>)) 
dt 

w (logy) 



-ew 2 (cos(a;t + (/>) + e(2sin 2 (o;t + . 



1)) 

(6) 

Plugging these approximations into (3) shows that u can 
be expressed approximatively as a polynomial of degree 
2 in esin(u;i + <p) and ecos(wt + <j)). Such polynomial 
computations can be easily done by a nonlinear analogue 
circuit. Figure 2 shows that the diode response remains 
almost the same when such polynomial approximations 
are used. Other simulations displayed on figure 3 show 
that such pre-compensation scheme are also robust to 
parameters uncertainties. 

5. CONCLUSION 

We have shown that it is possible to compute, via flatness- 
based motion planing techniques, adapted input current 
modulations in order to generate, in real-time and via 
analogue nonlinear circuits, a given modulation of the 
light entering the optic fiber. If the distortion along the 



Fig. 2. Input current modulations with e(t) given by (4): 
the solid curve corresponds to the response of the 
diode when u is perfectly related to y — y(l + 
e(t) cos(wt)) via (3); dashed line corresponds to the re- 
sponse of the diode when approximations (5) and (6) 
are used in (3). 

fiber is not negligible (long fiber and/or refractive index 
strongly dependent versus optic frequency), we have to 
consider also the phase dynamic of the light emitted by the 
diode (see, e.g., Agrawal [1997] for more details about the 
spectrum widening due to phase/amplitude variations). 
Following Lee et al. [2003] we have to complete (1) by 



dt 9 



aBr n I t h 



(X-l) 



where ip is the phase of the emitted light and a a pos- 
itive constant. Thus the emitted light entering the optic 
fiber is represented by the classical electric field (complex 
notations) 

E(t) oc \J P(t) exp(—vp(t)) exp(—iuj t) 

where u>q is the optical frequency (around 10 15 Hz). It is 
interesting to notice that the diode dynamics (described 
by three states (P, X, ip)) and one input I) is still flat with 
z = log P+FBt p t c P— — (p being the new flat output. This 

results from -^z = — 1+ B J pTc - . This structural property 
could certainly be exploited in other pre-compensation 
schemes taking into account such additional phase effect. 
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Fig. 3. Input current modulations with e(t) given by (4): 
as for figure 2, the solid curve corresponds to the 
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2 in order to test robustness versus large parametric 
errors. 
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